EFFECTS OF NON-CIRCULAR MOTIONS ON 
AZIMUTHAL COLOR GRADIENTS 



Eric E. Martfnez-Garcia 
Centro de Investigaciones de Astronomia, Apartado Postal 264, Merida 5101- A, Venezuela 

emartinez@cida.ve 
and 

Rosa A. Gonzalez-Lopezlira and Gilberto C. Gomez 

Centro de Radioastronomia y Astrofisica, UNAM, Campus Morelia, Michoacdn, Mexico, 

CP. 58089 

r . gonzalez@crya . unam . mx , g . gomez@crya . unam . mx 
ABSTRACT 



Assuming that density waves trigger star formation, and that young stars 
preserve the velocity components of the molecular gas where they are born, we 
analyze the effects that non-circular gas orbits have on color gradients across 
spiral arms. We try two approaches, one involving semi-analytical solutions for 
spiral shocks, and another with magnetohydrodynamic (MHD) numerical simu- 
lation data. We find that, if non-circular motions are ignored, the comparison 
between observed color gradients and stellar population synthesis models would 
in principle yield pattern speed values that are systematically too high for regions 
inside corotation, with the difference between the real and the measured pattern 
speeds increasing with decreasing radius. On the other hand, image processing 
and pixel averaging result in systematically lower measured spiral pattern speed 
values, regardless of the kinematics of stellar orbits. The net effect is that roughly 
the correct pattern speeds are recovered, although the trend of higher measured 
fl p at lower radii (as expected when non-cir cular motions exist but are ne glected) 
should still be observed. We examine the iMartmez- Garcia et al.l (120091 ) photo- 
metric data and confirm that this is indeed the case. The comparison of the size 
of the systematic pattern speed offset in the data with the predictions of the 
semi-analytical and MHD models corroborates that spirals are more likely to end 
at Outer Lindblad Resonance, as these authors had already found. 
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Subject headings: galaxies: kinematics and dynamics — galaxies: magneto- 
hydrodynamic simulations — galaxies: stellar content — galaxies: spirals - 
galaxies: structure 



Introduction. 



Until the detection of an azimuthal color gradient across one of the arms of the SAc 
galaxy M 99 (IGonzalez fc Grahamlll996l . GG96 hereafter), only sparse evidence of star for- 
mati on triggered by spiral density way es had bee n found in stel lar coun ts in the Milky 
Wav Jsitniklll989l . Il99ll : lAvedisovalll989l ) and M 3lJEfremovlll980allbl. Il985h. More recently- 



applyi ng the same method defined by GG96 and described below, iMartfnez- Garcia et al. 



(120091 . MG09 hereafter) examined a sample of 13 spiral galaxies of types A and AB, and 
found color gradients consistent with theoretical expectations in 10 of their objectso Al- 
though they did not compute stellar orbits, both GG96 and MG09 implicitly assumed that 
all stars, including those recently born near spiral arms, move in circular orbits. MG09 did 
investigate the effects of variable circular speeds and variable densities on observed color 
gradients and found them to be negligible. However, complicated non-circular motions have 
been reported in studie s about the migration of young stars following s t ar formation trig- 
gered by spiral shocks (lYuanlll969l ; IWielenlll979l ; Fernandez et al.ll2008l ). iBashl (Il98ll ) had 
also noticed that the ballistic orbits calculated from galactic HII region complexes follow 
non-circular trajectories that initially move along (and not across) the spiral arms. 

In the present work, we examine how the presence of non-circular motions would modify 
the pattern speeds derived from color gradients under the assumption of circular orbits. 

We will explore two different approache s. The first one (see § 31) involves the semi- 



anay litical solutions obtained for spiral shocks (jRobertslll969l ; IShu et al. 



1973: Gittins k Clarke 



20041 ): we assume that newly born stars preserve the orbital motion of the shocked gas where 
they form. The second approach (see § H]) is based on data from MHD simulations; we follow 
the gas flow vectors near spiral arms. We always assume that stars are triggered all along 
the studied regions of spiral arms in the respective model. Other methods, not discussed 
here, may involve orbit calculations for young stellar groups. 



1 See also 



Grosb0l fc Dottori 



(|2009f ) for an infrared study of young stellar complexes in NGC 2997. 
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2. The GG96 method. 



The photometric technique employed by GG96 and MG09 uses three optical bands, 
g, r and i, plus the near infrared J band (see tabled]). With these filters, the reddening 
insensitive and star-formation sensitive photometric index Q is obtained: 



Q{rJgi) = (r - J) - 0.82(g - i), (1) 

j2.05j2.50 

Q(r Jgi) = log 10 * i 05 - (2) 

r i 

In star-forming regions, the Q(rJgi) index has higher values, because the g and J bands 
in the numerator of equation [2] trace the light from red and blue supergiants, respectively. 
When the star- formation activity is poor, the value of the Q(rJgi) goes down. 

For each st udied galaxy and at each wavelength, the spiral arms are first unwrapped 



( live et al.l Il982l ) , by plotting them in a In R versus 9 map (see figures [fl [2]) . All the spiral 
arms in the disk galaxy are inspected, in search of color gradients. If a candidate region is 
found, the arm is straightened, by adding a different phase-shift to 9 at each value of In R, 
until the arm appears as a horizontal line (see figure [3]). This procedure allows us to add 
the data from different In R, and thus to increase the signal-to-noise ratio of the light profile 
as a function of 9. 1-D plots of Q vs. 9 are obtained from the unwrapped and straightened 
galaxy images. These are then compared to 1-D plots of (g — J), that trace the dust lane 
location, and to 1-D plots of the .fT-band (or K' or K s , depending on what we have available) 
data, that trace the density wave (under the assumption that the near-IR emission of red 
supergiants can be neglected). The regions that seem to match theoretical expectations 
according to these plots are selected and compared to stellar population synthesis models, 
that provide Q function of stellar age, t 



age- 



If one assumes that stars form in the site of the shock, and that they age as they move 
away from this location, then distance from the dust lane (at constant radius) parametrizes 
stellar age. In fact, stretching the model Q(t age ) to fit the data (where Q is a function of 

Table 1. Filters employed in the GG96 method 



Filter \ ef f FWHM 

g 5000A 830A 

r 6800A 1330A 

i 7800A 1420A 

J 1.25/im 0.29/xm 




Fig. 1. 



Schematic of a spiral disk. 
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Fig. 2. — Unwrapped version of figure [U Logarithmic spiral arms now appear as straight 
lines with slope = cot(-z), where i is the arm pitch angle. Horizontal lines are located every 
90°. Vertical lines correspond to radii marked with white circles in figure [H 
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Fig. 3. — Straightened version of figure [2j Now spiral arms appear as horizontal lines. The 
lni? coordinate remains the same, but the 9 coordinate has a phase-shift that depends on 
In R. 
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angular or linear distance) fixes the ratio between the distance from the shock, d, and the 
age of the stellar population. If, in addition, the rotational velocity is known, it is possible 
to find the angular velocity of the spiral pattern, Q p , as follows: 

fi P = Uot(R) - ■ (3) 

-^mean \ Cage / 

Here, R mea n is the mean radius of the region where the gradient has been detected, measured 
from the center of the galaxy; v rot (R) is the rotational velocity of the disk at such radius 
(obtained from the literature); d is distance from the shock, and £ age is the stellar age obtained 
from fitting the population synthesis model to the observations. With the pattern speed, we 
obtain the location of major resonances: the Outer Lindblad Resonance (OLR), corotation, 
and the 4/1 resonance. This resonance positions are compared with the observed spiral end 
points. 

As mentioned before, this procedure implicitly assumes circular motion for the stars 
involved, even though no stellar orbits have been computed. 



3. Semi-analytical shock solutions. 

Shock solutions o btained f rom s emi-analytical approaches have shown that gas stream- 
lines are not circular. Roberts! ( 1969 ). for example, found that streamlines appear as sharp- 



pointed ovals (see figure HJ). Although such solutions were found for the tightly w ound 



appro ximation, they provide a good reference for the study of real galaxies (see also ISaaf 



1974 ). Here, we obtain gas streamlines under the assumption that youn g stars follow the 



motio n of the molecular gas, and by adopting a procedure similar to that of lGittins fc Clarke 



□ Our adopted model has a pattern speed of 13 km s _1 kpc" 1 , and a flat rotation 
curve with v ro t = 220 km s _1 . From R = 5 kpc to R = 10 kpc, streamlines were calculated 
in steps of 0.05 kpc. Since the stellar models give us Q as a function of t age , we can assign 
a stellar age to each point in the streamlines; in order to do so, we need to assume that the 
onset of star formation occurs at a certain orbital time t, where t — corresponds to the 
shock position. 



2 See appendix 1X1 The reference system is the rotating frame of the spiral pattern. 
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1.8 2 2.2 

In [R(kpc)] 

Fig. 4. — Typical streamlines obtained from semi-analytical solutions of spiral, trailing type, 
shocks. Dotted line: spiral potential minimum; dotted circles: shock position. Average radii, 
from left to right, are 6, 7, 8, 9, and 10 kpc. 
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3.1. Star formation onset delay. 



Martmez-Garcfa et al.l (120091 ) use the term "timescale for star formation" as the duration 
of the burst, that is estimated to be 2 x 10 7 years. They do not make any assumptions 
about the time it takes for star formation to begin after the shock. A preliminary inspection 
of their findings shows that the shock (as traced by the dust lanes) and the onset of star 
formation (marked by t age = in the stellar models) are not always at the same location. 
In real galaxies, whether there is a delay between the shock and the star formation onset 
is not well understood yet. In fact, the time needed for a diffuse cloud of neutral gas to 
first become a dense cloud, then a molecular cloud, and finally a self-gravitating cloud is 
10 7 years. Once this cloud is formed, the onset of st a r formation may b e very fast 



(jVazquez-Semadeni et al.l 120071 ; iHeitsch fc Hartmannl 120081 ) . lEgusa et al.l (120091 ) , based on 
an observational study between the peaks of molecular (CO) and young stellar (H a ) arms, 
report a delay f o r the onset of star formation of ~ 5 — 30 Myr in 5 galaxies (out of 13); 
Tamburro et al.l (120081 ) estimate a timescale for star formation in the range 1-4 Myr, from 
a study of angular offsets between HI and 24/im emissivity peaks in a sample of 14 disk 
galaxies % 

In order to test for the effects of a star formation onset delay on azimuthal color gra- 
dients, we try three cases. The first one has a delay (i.e., the time it takes for the gas to 
move from the shock to the location where star formation begins) St s h oc k — yr; the other 
two have St s h oc k = 1 x 10 7 yr and 5t s h oc k = 2 x 10 7 yr, respectively. 



4. MHD simulations. 



The gas simulations were performed with a version of the ZEUS code (jStone fc Norman 



1992al lbL that is an eulerian, time-explicit, finite-difference code for ideal MHD simulations. 



We employ a 2D grid with 500 x 500 points in polar coordinates. The radial extent goes 
from R ~ 3 to 30 kpc, and the azimuthal one from 9 ~ to tc radians. We assume that 
the half disk simulation data have a 180° "mirror" symmetry. The gaseous disk follows 
an isothermal equation of state with a temperature of 10900° K. At the beginning of the 
simulation, the magnetic field has a toroidal geometry and a value of 5 fiG at 8.5 kpc from the 
galactic center, although it rapidly evolves away from this setup. No self-gravity is included. 

The background gravitational potential has two components: one axisymmetric, that 
consists of a bulge, a disk, and a halo, while the other is the non-axisymmetric spiral arm 



3 For these authors, the term "star formation timescale" is equivalent to star formation onset delay. 
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perturbation. The adopted bul ge and disk are described in I Allen fc Santillanl (Il99ll ). The 
dark matter halo is NFW type (INavarro et al.lll996l . 119971 ). with density: 



p{R) 



Pn 



R/a h (1 + R/a h ) 



2 ■ 



(4) 



where p & =1.021 x 10~ 2 M & /pc 3 and a fc =15.133 kpc. The two -arm spiral potential has a 
pitch angle of 15.5°, and is self-consistent in the stellar orbits sense (IPichardo et al.ll2003l ). 
The simulation is performed in this stellar arms reference frame, that rotates with velocity 
a p = 20 km s- 1 kpc' 1 . 

We allow the simulation to evolve during 8 x 10 8 yr (see figure [5]). By this time, the gas 
develops spiral arms as a response to the imposed perturbation, initial transients subside, 
and changes in the simulation are only observable at very long time scales. Therefore, we 
assume that density and velocity distributions at this time represent a close approximation 
to a steady state solution. 

As shown in figure the gas in the simulat ion has responded to the spiral perturba- 
tion with a four arm pattern (IMartos et all 120041 ) . A hypothetical circular gas orbit would 
encounter four shocks before completing its transit around the disk. In this case, the cor- 
responding velocity field would be very different from the case where only two shocks are 
considered. For the present investigation we focus on the main shock (i.e., the one closer 
to the potential minimum of the spiral perturbation), although spiral-arm triggering of star 
formation may also take place in the secondary shocks. 

In order to obtain the orbit a parcel of gas would follow, we locate the main shock 
position by searching for the gas density maxima closest to the (stellar) spiral arms. These 
are taken as the start point for the integration of the gaseous orbits. The orbits are then 
calculated from the velocity data (in the spiral pattern reference frame), using a Runge- 
Kutta method, every AR = 0.005 kpc, and At = 7 x 10 12 s. A sample of orbits are shown in 
figure El Other gas orbits wi t h MH D simulations using the full time evolution of the gas can 
be seen in lGomez fc Martosl (120091 ). who find a close resemblance between gas orbits and the 
central family of ballistic orbits in stellar dynamics. 
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Fig. 5. — MHD simulation at 800 Myr. Colors represent the gas density, while arrows 
indicate the velocity field (in the rotating frame of the spiral pattern). The arrow in the 
lower right corner corresponds to ~ 80 km s . Corotation is located at R = 10.9 kpc. 
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Fig. 6. — Gas orbits for the simulation data, in the lni? vs. © plane. Initial radii (from left to 
right) are 5.5, 6, 7, 8, and 9 kpc. Dotted line: stellar arms density maximum; short-dashed 
line: primary gas arms density maximum; long-dashed line: secondary gas arms density 
maximum. 
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5. Analysis and Results. 



5.1. Q profiles. 

From the first approach, involving semi-analytical solutions, we obtain Q index profiles 
for the non-circular and circular cases, respectively, calculated in a narrow annulus with 
mean radius ~ 7 kpc; we show them in figure [TJ To produce the Q profiles, we assume a 
background population of old stars with an age of 5 x 10 9 yr (i.e., the density wave), plus 
a young burst of star formation lasting 2 x 10 7 yr. Both the old and the young populations 
have a Salpeter IMF, with Mi ower = O.1M and M upper = 10M Q @ For the present case, 
we adopt a constant fraction of young stars of 2% b y mass (see equation 9 in MG09)H We 
then use the stellar population synthesis models by Chariot fc Bruzual ( 20071 ) to find the 
changing g, r, i, and J emission of the population mixture as it evolves. 

In the circular case, the stellar complexes that give origin to the Q profile were all born 
at the same (\nR,9) position. In the non-circular case, due to the trajectories of the gas 
streamlines, the stars that give rise to the Q profiles were not all formed at the same location. 
The age of the stellar complexes that contribute to the Q profile as a function of distance 
from the shock position is plotted in figure [SJ As expected, young objects tend to spend 
more time concentrated toward the spiral shock (d = kpc), with the consequence that the 
peak of the Q profile (dashed line in figure [7J occurs closer to the shock. 

Another aspect to be noticed in figure [7J is that there is no "downstream fall" of the 
gradients (i.e., Q does not fall below, in this case, ~ 1.56 , contrariwise to what is sometimes 



observed with the data). iMartinez-Garcia et al.l (120091 ) had hypothesized that the fall of 
the observed Q profiles below the model values (assuming pure circular orbits) might be 
caused by stellar non-circular motions in the data. How ever, results may differ if ba llistic 



trajectories with postshock velocities are considered (e.g. iBash et al.lll977l ; lBashlll979l ) 



4 There is a conspicuous inverse correlation between the detection of azimuthal color gradients and the 
presence of HII regions (see MG09). 

5 The effects of variable densities of both the young and old populations on the color gradients have already 
been discussed bv lMartmez-Garcia et al.1 (|2009l ). Their main conclusion is that variable stellar densities can 
produced deformations in the expected color gradients. However, the estimated error introduced by these 
deformations is lower than the computed random error contributed by the combined uncertainties in the 
inclination angle, the rotation velocity, and the distance to the galaxy. 



-14- 



5.2. Pattern speeds. 

From the Q profiles, it is possible to derive pattern speeds by comparing stellar popu- 
lation synthesis models (Q vs. t age ) to the "data" (Q vs. d), assuming implicitly that stars 
move in circular orbits, as was done by GG96 and MG09, and was sketched above in section 

m 

In the case of the semi-analytical approach, "data" Q profiles were obtained from narrow 
annular regions, ~ 0.001 kpc wide. Given that we know the input pattern speed, we can 
search for systematic effects in the determination of Q p . Figure [9] shows the input pattern 
speed (long- dashed line), and Q p values that are obtained from the color gradients for three 
different presuppositions about the delay of star formation onset after the shock (0, 10, and 
20 Myr). For reasons explained below in § [6], rather than stretching the stellar population 
synthesis model to the "data", we just compare the positions of the maxima of the model 
Q, on the one hand, and of the shock location, on the other. 

There is a significant systematic effect, whereby the derived Q p is always larger than 
the input pattern speed; moreover, the effect decreases with galactocentric radius. Also, 
the difference between the input and the output pattern speeds decreases for larger star 
formation onset delays. 

In the case of the MHD simulations, we again analyze the pattern speeds derived when 
assuming a circular motion dynamic model. In order to obtain fl p , we compare the peak value 
of Q in the stellar population model to the density maximum of the gas in the simulations. 
In figure [TUl we show the fl p values obtained for 5t s h oc k of 0, 10, and 20 Myr. 

6. Pixel averaging 

One of the caveats of the GG96 method, that was not discussed in MG09, comes from 
image processing. Each pixel in the optical and infrared images actually contains information 
about many young stellar orbits that all fall within the same spatial region. The process 
of "unwrapping" the spiral arms then averages pixels in the 9 and lni? directions (the 
"straightening" of the arms is just a shift in the 9 direction and does not involve pixel 
averaging). In order to increase the signal-to-noise ratio of the Q profiles, an additional 
averaging in the lni? direction is done in selected regions. This last step averages together 
orbits with different angular speeds, and hence stars of slightly different ages at a fixed 
distance from the shock. 

Figure [TT1 shows the behavior of Q profiles from an unwrapped and straightened synthetic 
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Fig. 7. — Q index profiles inside a narrow annulus with mean radius R = 7 kpc, vs. azimuthal 
distance relative to the shock (d—0). Solid line: circular case; dashed line: non-circular case. 
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Fig. 8. — Stellar population ages inside a narrow annulus with mean radius R = 7 kpc, vs. 
azimuthal distance from the shock (d—0). Solid line: circular case; dashed line: non-circular 
case. 
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Fig. 9. — fi p values obtained at different galactocentric radii, from the numerical semi- 
analytical solutions of spiral shocks, under the (false) assumption that stars move in circular 
orbits. Solid line: star formation onset delay 5t s h oc k — yr; dotted line: St s h oc k — 10 Myr; 
short-dashed line: St s h oc k = 20 Myr. Long-dashed line: spiral perturbation input angular 
velocity. 
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R (kpc) 

Fig. 10. — Q p values obtained from the MHD simulation data, at different galactocentric 
radii, assuming a circular motion dynamic model. The "data" Q profiles were calculated in 
narrow annuli 0.001 kpc wide. Symbols as in figure [9j 
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image (see below), under the assumption of pure circular motion; the different profiles are 
taken from various mean radii, after averaging in InR. The shift of the peak toward d = is 
expected as we approach the corotation radius. However, due to the pixel averaging, a drop 
in the maximum value of Q is also obtained. If, in order to match the observed or synthetic 
profile peak, a constant downward shift is applied to the stellar population synthesis model, § 
the "wings" of the profiles will not match, since those of the stellar model will now lie within 
the "wings" of the data. This mimics wider gradients, and in order to fit the observed or 
synthetic profile we will need to overstretch the model. Consequently, the pattern speed 
obtained via equation [3] will be lower, because the stellar age (t age ) at a fixed distance d 
from the shock will be underestimated. The effect amounts to ~ 1 to 2 km s~ x kpc -1 for 
moderate pixel averaging, but can be as large as 10 km s~ x kpc -1 when the averaged region 
in In R is very extended. 

One way out of this problem is to compare only the distances between the onset of 
star formation and the Q peak, in the model and the data, instead of the whole Q profiles, 
since image processing affects less the peak positions than the profile shapes. We have used 
this approach with our semi-analytical calculations and MHD simulations; unfortunately, 
in the case of real data, the position of the star formation onset is generally unknown (see 
also § EZEp . 

To better estimate the effects of averaging, with the R, 9 and t (time elapsed since the 
shock, t=0 corresponds to the shock location) data of the streamlines obtained from the semi- 
analytical solutions, we constructed images of a synthetic "galaxy" as it would be observed in 



the q , r, % and J bands with the aid of stellar population synthesis models (ICharlot fc Bruzual 



20071 ) (the fraction of young stars constitute is 2% by mass). The model face-on galaxy is 
at a distance of 35 Mpc and the synthetic images have 1" x 1" pixels. Whenever various 
streamlines were located at the same pixel, their elapsed times since the shock were averaged. 
Next, the synthetic images were unwrapped. Each hyperpixel of the unwrapped images is 
3° long in the 9 axis; the InR axis has 250 hyperpixels in total, from the center of the 
galaxy to R — 10 kpc. Then, the synthetic unwrapped images were straightened according 
to the shock's pitch angle (~ 5.266°), that is smaller than the potential minimum's pitch 
angle (~ 5.739°; see appendix EJ. Finally, Q profiles were obtained from these "images" in 
sections with a width of 5 hypercolumns in In R. Pattern speeds were measured by comparing 
the "observed" Q (vs. d) from the simulated images, at different radii, with Q (vs. t age ) from 



6 Vertical shifts are also expected if the metallicities, or the ratios of young to old stars of model and 
data are different; these differences do not significantly affect the derived pattern speed (MG09), because 
the horizontal extension of the profiles is not noticeably affected. 



-20- 



5 10 15 




5 10 15 

d (kpc) 



Fig. 11. — Q profiles obtained from unwrapped and straightened images at different mean 
radii, after averaging in lni?. The dynamic model assumes pure circular motions with 
v ~ 220 km s" 1 , and Q p = 13 km s _1 kpc -1 . -Rmean denotes the region mean radius, d is 
the azimuthal distance fr om the shock {d — 0). Stellar population synthesis models from 
Chariot fc Bruzuall ((2003), with IMF lower limit of 0.1 M , and upper limit of 10 M . 



stellar population synthesis models 



In figure [12j we illustrate the effects of averaging. The dotted line is the input pattern 
speed. If the synthetic galaxy has circular orbits and the galaxy is analyzed assuming 
circular orbits, then one obtains the fl p values depicted by squares. The solid black squares 
are obtained with the GG96 "stretching" method, whereas the comparison between model 
and data Q maxima positions relative to the shock yields the open symbols. Clearly, the 
solid black squares are biased while the open ones are not. 

On the other hand, if the synthetic galaxy has non-circular orbits but is analyzed under 
the assumption of circular orbits, then one obtains the Q p values shown as triangles. The 
solid black triangles are the values obtained with the GG96 method; the empty triangles are 
the pattern speeds found when comparing only the positions, relative to the shock location, 
of model and "observed" Q peaks; and the solid line represents the biased values measured 
from the semi-analytical solutions, without averaging, if existing non-circular motions are 
ignored. The reason for this bias towards higher measured pattern speeds is that real gas 
streamlines in a steady rotating spiral shock turn somewhat along the arm after passing 
through the shock. Consequently, stars close to the shock are slightly older than would 
be expected in a circular model with the same fl p . The observer using a circular model 
would infer that the gas flow into the arms is smaller than is actually the case, i.e., that 
the difference between the stellar orbital velocity and the pattern speed is smaller than in 
reality; inward of corotation, this means that the observer would overestimate the pattern 
speed. Also inside corotation, the effect decreases with galactocentric radius, as the shock 
strength diminishes. 

The open triangles product of the peak comparison are much less sensitive to the compet- 
ing systematics introduced by orbit averaging. However, ironically, the original "stretching" 
method (solid black triangles) seems much less biased towards higher values if non-circular 
motions are neglected, although the dependence of the bias on galactocentric radius is still 
noticeable. This is a fortunate turn of events, because (1) in real images it is hard to pinpoint 
the location of the star formation onset, and (2) the radial dependence of the bias allows us 
to detect the presence of non-circular motions and to confirm the link between star formation 
and disk dynamics. 



7 The adopted models throughout this investigation have an IMF upper limit of 10 Mq. Models with 
higher upper limits would peak at younger ages, closer to the shock; see MG09 for models with an IMF 
upper mass limit of 100 Mq. 
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Fig. 12. — Q p measurements. All the analysis methods assume circular kinematics, re- 
gardless of the actual kinematics of the input galaxy. Dotted line: input pattern speed. 
Squares: measurements, with pixel averaging, for a synthetic galaxy with purely circular 
motions obtained through the GG96 "stretching" method (solid) and via the comparison 
of Q peaks (open). Solid line: Q p values derived without averaging in the case where the 
synthetic galaxy has streamline trajectories described by the semi-analytic model. Trian- 
gles: measurements, with pixel averaging, for a synthetic galaxy with streamline trajectories 
described by the semi-analytic model deduced from the GG96 method (solid) and with the 
Q peak method (open). 
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Discussion. 



From theory (e.g.. lGittins fe Clarkdl2004l ). it is well known that the spiral shock strength 
diminishes as one approaches corotation. The shock weakening implies fewer radial move- 
ments and more circular trajectories. Figures [9] (corotation radius, Rcr, located at ~ 17 
kpc) and [TO] (Rcr ~ H kpc) support this conclusion; they also allow us to predict that, in 
general, Q p values derived from color gradients observed inside corotation in real data will 
be too high if a circular dynamic model is adopted^ For color gradients observed beyond 
corotation, we expect again to overestimate the pattern speed, but the difference between 
the actual and the measured Q p will now grow with radius. In this case, the radial velocities 
of the streamlin es after the sh ock have an outward, rather than an inward, component (see 
leading case in lRobertslll969l ). If, however, spiral shocks are not as strong as inside corota- 
tion, inverse post-corotation gradients will yield unbiased pattern speed measurements. 

How would these results affect the Martmez-Garcia et al. J2009h conclusions about the 
end points of spiral patterns? Their figure 4 displays the ratio, -R res /-R e nd, of the resonance 
radii (calculated from their pattern speed measurements and rotation curves) to the observed 
spiral end points in their deprojected, near-infrared galaxy images. If their Q p measurements 
(most of them inside corotation) are systematically too high, then the resonance radii, and 
their ratio to the spiral end points would be systematically smaller than the true values. 

In order to better quantify the impact of non-circular motions on the MG09 results, 
we define the expression 5Q P = (fidata — ^ P )/^ P - Here, fidata is the pattern speed obtained 
from the data, by comparing the observed color gradient candidates with the stellar models, 
under the assumption that stars move in purely circular orbits; Q' is the pattern speed that 
the spiral should have if it ends at certain resonance. For the 4:1 resonance, for example, we 
have: 



p 



Vrot 



R 



1 - 



y/2 



(5) 



1:1 



where v TO t, is the circular orbital velocity (obtained from the literature, see MG09). And if 
the spiral ends at the OLR: 



8 For some radii (e.g., ~ 6 and 8.5 kpc) in the MHD simulation, the measured tt p will be actually slightly 
lower (~ 1 km s -1 kpc -1 ) than the input pattern speed of 20 km s -1 kpc -1 . 
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To calculate 5Q p for the 23 regions analyzed in MG09, we take both fi^ata and -Rend from 
their Table 4. In figures [TBI and [T4"l we plot 5Q P vs. -R m ean/-Rend, where -R me an is the mean 
orbital radius of the studied region. Figure [13] co rresponds to the case wh e re the end point 



of the spiral pattern is fixed at the 4:1 resonance. IContopoulos fc Grosboll (119861 ) state that 
strong spirals (Hubble types Sb or Sc) may be truncated at this resonance. The 6fl p values 
show the trend with radius expected from figures M and [TDJ This would imply that we are 
actually detecting non-circular motions from the photometric data! 

However, three discrepancies with respect to theoretical expectations are also present. 
The first one is related to the systematic effects discussed in appendix EJ In spite of the 
tendency to actually measure values around the true pattern speed with the method used 
by GG96 and MG09 (solid black triangles in figure [T2l . there is a lack of points below the 
Sflp = line. The second difference is the gap around i? mean /i? e nd ~ 0.6, where a point 
with huge vertical error bars is located^ If spirals are really truncated at the 4:1 resonance, 
and hence inside or near corotation, one would expect a continuous distribution of points 
between the ILR and the spiral end point (i.e., -R mea n/-Rend = 1)- The third inconsistency 
concerns the actual values attained by 5Q P at low -R mG an/-Rend- From the solid black triangles 
in figure [12], the maximum expected value of SQ P is < 0.5. The high 5Q P values would imply 
the presence of extremely strong shocks near the spiral arms. 

Figure [T4"] in contrast, shows the behavior of SQ p for the case when the spiral end points 
are fixed at the OLR. If spirals stop near this resonance, the expected -R m ean/-Rend value for 
corotation is ~ 0.59, and the gap in the data distribution is naturally explained. We also 
observe now that some points are located below the 5Q P = line, as expected from figure [12] 
owing to the systematic error introduced by the "stretching" GG96 method and discussed 
in § El The three last points, marked with a long-dashed ellipse, represent the "inverse color 
gradients" (i.e., located beyond corotation) examined in MG09 (these points were treated as 
being inside corotation in figure [13]). The points surrounded by a short-dashed ellipse belong 
to the galaxy NGC 578. If NGC 578 actually ends before or at corotation (as argued by 
MG09), its position in figure [TBI is in concordance with theoretical expectations (i.e., close to 
the Sflp = line). On the other hand, NGC 578 may be the only gal axy in the MG09 sample 



with a spiral pattern that stops at the 4:1 resonance, as predicted by IContopoulos &: Grosbol 



9 MG09 argued that this feature could be due to a star formation event, not triggered by spiral waves, 
located near corotation in NGC 1421. 
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8. Conclusions. 

Under the assumption that the orbits of young stars preserve the velocity components 
of the parent molecular clouds where they form, we have analyzed the effects that non- 
circular motions would have on azimuthal color gradients. Semi-analytical calculations and 
MHD simulations show that the spiral pattern speeds derived from the comparison between 
color gradients and stellar population synthesis models, assuming purely circular motions, 
would have values systematically higher than the real ones for regions within corotation. 
Also inside corotation, the effect would decrease with galactocentric radius. Non-circular 
motions, however, would not prevent the detection of legitimate azimuthal color gradients 
in real galaxies. 

Using a synthetic image, we have also analyzed the effects of image processing and pixel 
averaging on the method applied in GG06 and MG09 to detect color gradients and derive 
pattern speeds. We have found that pixel averaging (due to image processing) systematically 
decreases the derived Q p , such that it nearly compensates for the systematic effect introduced 
when neglecting non-circular motions in the analysis. The net result is that the correct spiral 
pattern speeds and resonance locations can be obtained. Nevertheless, a residual trend of 
slightly higher pattern speeds at lower radii can still be discerned (solid triangles in Fig. [12] 
and Fig. [T3j), so that it is possible to detect the presence of non-circular motions and confirm 
the link between star formation and disk dynamics. 

We have re-examined the results obtained by MG09. When normalizing the mean radii 
where gradients were found by the end point radii, in order to treat the whole sample as a 
single galaxy, we were able to reproduce the trend of 5Q P with radius expected if non-circular 
motions are ignored (as these authors did). The size of the observed 5Q p in the data only 
matches the theoretical expectations if spiral patterns end at the OLR, and not at the 4:1 
resonance (cf. Fig. [12] with both Fig. [T31 and Fig. |T4"|) . 

Future studies of azimuthal col or gradients in disk galaxies by orbit calculations should 



include the effects of disk heating (lAsiain et al.lll999l ). and a correct modeling of the IMF 



as a discrete distribution. These factors may be important, since high mass stars may take 
different trajectories when co mpared to low mass stars, due to the effects of dynamical 



friction (jChandrasekharl Il943l ) . 



We thank the anonymous referee for his/her very constructive suggestions and comments 
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Fig. 13. — 5Q P (see text) vs. -R m can/-Rend, for the 23 regions in Table 4 of MG09. The spirals 
are assumed to end at the 4:1 resonance. The dotted line indicates 5Q P = (i.e., no difference 
between the real and the measured pattern speeds). 
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Fig. 14. — Same as figure [13] for the case where spirals are assumed to end at the OLR. 
Long-dashed ellipse: Surrounds points that presumably correspond to inverse color gradients, 
outside corotation. Short-dashed ellipse: Surrounds points from regions in NGC 578, whose 
spiral pattern might end at or inside corotation (MG09). 
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A. Streamline locus determination from semi-analytical solutions. 



The differential equations needed to obtain the density and ve l ocity c omponents, prod - 
uct of spiral shocks, ar e fully explained in iRobertsI (119691 ) , IShu et al.l (119721 ) , IShu et al.l (119731 ) , 
and iGittins fc Clarkel ( 2004T ) J 10 l The independent parameters used in this investigation for 
the solutions are: number of arms, m = 2; spiral arms pitch angle % such that sini = 0.1; 
angular speed of the spiral pattern, Q p = 13 km s _1 kpc -1 ; effective speed of sound, a = 8 
km s _1 ; amplitude ratio of the perturbed spiral field to the axisymmetric field, F = 0.05; 
rotation velocity of materialS v rot = 220 km s l . In figure [151 we show the positions of the 
sonic point and the spiral shock relative to th e potential m ini mum (n = 0), with the pa ram- 



eters described above. As already noticed by iTosal (119731 ) and lGittins fc Clarkel (120041 ) . the 
shock front moves away from the potential minimum as the radius increases. In this case, 
for R ~ 5 kpc, the shock front is slightly ahead of the potential minimum. These relative 
positions may change for another set of parameters. 

The solutions are obtained in curvilinear coordinates r] (perpendicular to the spiral 
equipotential curves) a nd $ (para ll el to the s piral equipotentia l curve s) . We adopt the rj 



and £ definitions from IShu et al.l (11973 



and IGittins fe Clarkel (12004 ) . These definitions 



differ from those used in IRobertsI (119691 ) and IShu et al. <jl972h by a multiplicative factor 
mj sini, and by an additive constant, chosen such that rj = determines the location of 
the potential minumum. The corresponding expressions in terms of the radius, R, and the 
angular coordinate in the rotating frame of the spiral pattern, O = 9 — Q p t, are as follows: 



m 



tan % \R a 



+ m(0 - tt p t) + 7T, 



(Al) 



10 See also lFuiimotol (|1968l ). and llshibashi & Yoshiil (|l984l) . 

11 A generic flat rotation curve is adopted, such that the epicyclic frequency K = We also assume 

that F ^constant for the selected range of radii (5-10 kpc). 
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e = _ mhl (^ + J^ {6 -n p t). (A2) 

\rt a J tan 2 

The rj coordinate is obtained from the shock solution, together with the velocity com- 
ponents = u^o + Uq, and u v = u v0 + u n \. The subscript labels the unperturbed velocity 
(i.e., in absence of a spiral perturbation), and the subscript 1, the perturbation due to the 
spiral gravitational field. The £ coordinate is obtained by solving the equation: 



that can also be expressed as: 



f = (A3) 
dr] u v 



— ) dr] = ^stream ~ £sub- (A4) 



Vsub 



dr] 



The shock on the subsonic branch occurs at r] sub and £ su &. r] stream and Stream are coordinate 
values at any point along the streamline; their maximum possible values are r] sup and £ snp , 
i.e., the point at which the shock occurs on the supersonic branchial 



In order to obtain £ sub , we proceed in the way described below. We define the quantity 

x = R 



r e 'su P ( _r V ,0/ 
B' — B' 

y sup ^ sub J 



(A5) 



where an d @' are the (7^) an d @ values obtained when £ sub = and _R a = 1. i? is 

the average radius of the streamline. 



We use the x value to get: 

/ Tfl \ 

isub = - m + - — 2- In x. (A6) 
V tan 1 J 

With this definition, the mean value of (^j^j w iU De Rq- 

The time elapsed since the shock, t, is obtained by integrating the velocity along the B 
direction: 



Vsup — Vsub = 27T. In theory, £ sup — £ su h = but, du e to a slight nonclosure of the streamlines with 
radius, the numerical values may differ fsee lshu et aDll97il . 
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[ @ °~ RdO 

t= — , (A7) 

Je sub u e 

where t = when 6 stream = 6 snfe , and 



tt„ + 1i£ cot 2 , . r 

Me = ! . (Ac 

cot « cos i + sin i 



REFERENCES 

Allen, C, & Santillan, A. 1991, Revista Mexicana de Astronomia y Astrofisica, 22, 255 

Asiain, R., Figueras, F., & Torra, J. 1999, A&A, 350, 434 

Avedisova, V. S. 1989, Astrophysics, 30, 83 

Bash, F. N., Green, E., & Peters, W. L., Ill 1977, ApJ, 217, 464 

Bash, F. N. 1979, ApJ, 233, 524 

Bash, F. N. 1981, ApJ, 250, 551 

Chandrasekhar, S. 1943, ApJ, 97, 255 

Chariot, S., Bruzual, G. 2007, private communication 

Contopoulos, G., & Grosbol, P. 1986, A&A, 155, 11 

Efremov, Y. N. 1980, Soviet Astronomy Letters, 6, 152 

Efremov, Y. N. 1980, Soviet Astronomy Letters, 6, 184 

Efremov, Y. N. 1985, Soviet Astron. Lett., 11, 69 

Egusa, F., Kohno, K., Sofue, Y., Nakanishi, H., & Komugi, S. 2009, ApJ, 697, 1870 

Fernandez, D., Figueras, F., Torra, J. 2008, A&A, 480, 735 

Fujimoto, M. 1968, IAU Symposium, 29, 453 

Gittins, D. M., & Clarke, C. J. 2004, MNRAS, 349, 909 



-31 - 




Fig. 15. — Sonic point and shock positions at different radii. Dashed line: rj sp , sonic point 
position relative to the potential minimum (dotted line, rj = 0); solid line: rj su b, position of 
the spiral shock relative to the potential minimum. 



-32 - 

Gomez, G. C, & Martos, M. A. 2009, Chaos in Astronomy Conference 2007, Astrophysics 
and Space Science Proceedings, Edited by G. Contopoulos and P.A. Patsis. 

Gonzalez, R. A., & Graham, J. R. 1996, ApJ, 460, 651 (GG96) 

Grosb0l, P., & Dottori, H. 2009, A&A, 499, L21 

Heitsch, F., & Hartmann, L. 2008, ApJ, 689, 290 

Ishibashi, S., & Yoshii, Y. 1984, PASJ, 36, 41 

lye, M., Okamura, S., Hamabe, M., & Watanabe, M. 1982, ApJ, 256, 103 

Martinez- Garcia, E. E., Gonzalez-Lopezlira, R. A., & Bruzual-A, G. 2009, ApJ, 694, 512 

Martos, M., Hernandez, X., Yanez, M., Moreno, E., & Pichardo, B. 2004, MNRAS, 350, L47 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Pichardo, B., Martos, M., Moreno, E., & Espresate, J. 2003, ApJ, 582, 230 

Roberts, W. W. 1969, ApJ, 158, 123 

Shu, F. H., Milione, V., Gebel, W., Yuan, C, Goldsmith, D. W., & Roberts, W. W. 1972, 
ApJ, 173, 557 

Shu, F. H., Milione, V., & Roberts, W. W., Jr. 1973, ApJ, 183, 819 

Saaf, A. F. 1974, ApJ, 189, 33 

Sitnik, T. G. 1989, Soviet Astron. Lett., 15, 388 

Sitnik, T. G. 1991, Soviet Astronomy Letters, 17, 61 

Stone, J. M., & Norman, M. L. 1992, ApJS, 80, 753 

Stone, J. M., & Norman, M. L. 1992, ApJS, 80, 791 

Tamburro, D., Rix, H.-W., Walter, F., Brinks, E., de Blok, W. J. G., Kennicutt, R. C, & 
MacLow, M.-M. 2008, AJ, 136, 2872 

Tosa, M. 1973, PASJ, 25, 191 



-33- 



Vazquez-Semadeni, E., Gomez, G. C, Jappsen, A. K., Ballesteros-Paredes, J., Gonzalez, 
R. F., & Klessen, R. S. 2007, ApJ, 657, 870 

Wielen, R. 1979, The Large-Scale Characteristics of the Galaxy, 84, 133 

Yuan, C. 1969, ApJ, 158, 889 



This preprint was prepared with the AAS IATgX macros v5.2. 



